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Abstract 

Spherically symmetric static solutions of the Einstein equations 
with a positive cosmological constant for the energy-momentum tensor 
of a barotropic perfect fluid are governed by the Tolman-Oppenheimer- 
Volkoff-de Sitter equation. Sufficient conditions for existence of solu¬ 
tions with finite radii are given. The interior metric of the solution 
is connected with the Schwarzschild-de Sitter metric on the exterior 
vacuum region. The analytic property of the solutions at the vacuum 
boundary is investigated. 


1 Introduction 


We consider a static and spherically symmetric metric 

ds 2 = g ilv dx il dx v = e 2F{r) c 2 dt 2 - e 2H ^dr 2 - r 2 (d6 2 + sin 2 Odcf 2 ) 
which satisfies the Einstein-de Sitter equations 

7? 1 p a _8 ttG 
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for the energy-momentum tensor of a perfect fluid 

T ]W = ( c 2 p + P)U^U^ - Pg» v . 


Here R^ u is the Ricci tensor associated with the metric g^dx^dx 1 ', R = 
g a ?R a p is the scalar curvature, and c, G are positive constants, the speed 
of light, the gravitational constant. A is the cosmological constant which is 
supposed to be positive, p is the mass density, P is the pressure and £/ M is 
the four-velocity. See [TU] §111]. 

Historically speaking, the cosmological constant A of the above Einstein- 
de Sitter equations was introduced by A. Einstein [7], 1917, and was discussed 
soon by W. de Sitter [16j. Although it was introduced for a static universe, 
it was not necessary for an expanding universe. So, later Einstein wrote to 
H. Weyl “If there is no quasi-static world, then away with the cosmological 
term! ” on May 23, 1923, and finally he rejected the cosmological term in [8], 
1931, for the reason that it is not necessary to explain the Hubble’s report 
on the redshifts of galaxies showing the expansion of the universe. (See [BEl 
§15e].) However, although the original motivation for introducing the cos¬ 
mological term disappeared, its status in the cosmological theories remained 
and it revived with new meanings in the recent development of the theories 
and observations. For the details, see the review [3J 

For the static spherically symmetric metric, the Einstein-de Sitter equa¬ 
tions are reduced to 



The coefficients of the metric are given by 


g2 F(r) _ k+£ - 2u(r)/c 2 ^ 

u(r), n + being the function and the constant specified later, and 



2 Gm(r) A 2 


—r . 
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When A = 0, the equations (jT() turns out to be 


dm 
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and this ([2]) is called the Tolman-Oppenheimer-Volkoff equation. It was 
derived in [13], 1939. Therefore we shall call (HD with A > 0 the Tolman- 
Oppenheimer-Volkoff-de Sitter equation. In this article we investigate 
this Tolman-Oppenhcimer-Volkoff-de Sitter equation (]T|) . Throughout this 
article we keep the following 

Assumption The pressure P is a given function of the density p > 0 
such that 0 < P and 0 < dP/dp < c 2 for p > 0 and P —> 0 as p — * +0. 
Moreover we assume that there are positive constants A, 7 and an analytic 
function Q on a neighborhood of [ 0 , + 00 ) such that fi( 0 ) = 1 and 

P = Ap' y Q(Ap 1 ~ 1 /c 2 ). 

We assume that 1 < 7 < 2. 


Actually we are keeping in mind the equation of state for neutron stars: 


P = Kc 5 


ft q 4 dq 

L \A + q 2 ' 


p = 3 Kc 3 / \/l + q 2 q 2 dq , 


K being a positive constant. See p2], [IT, p. 188]. 


The Tolman-Oppenheimer-Volkoff-de Sitter equation has already been 
systematically investigated from the physical point of view by C. G. Bohmer, 
n. Pi- However in the study by C. G. Bohmer there is supposed to exist a 
positive density pb, so called the ‘boundary density’, at which the pressure P 
vanishes, that is, P > 0 yy p > pb, and cosmological constants A satisfying 

A < —— pb are considered. This situation is not treated in this article. 

<r 


(3) 


For the sake of notational conventions, we shall denote 

2 Gm A 


K(r,m ) := 1 


c 2 r 


—r 2 , 

3 ’ 


and 

(4) Q(r,m,P) := G(m + ^-P^j -^r 3 , 

so that the second equation of (JT[) reads 


d r = - {p + P/c ‘ ) 9pp. ih). 

dr r z K{r , m) 
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We consider the equation (jT]) on the domain 


D = {(r, m, P)|0 < r, \m\ < +oo, 0 < p, 0 < n{r, m)}. 


Now we have a solution germ (m(r),P(r)) at r = +0, given the central 
density p c > 0 with P c := P(p c ), such that 


(5) m = —p c r 3 + 0(r 5 ), 
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as r —> +0. Proof is the same as that of [121 Proposition 1], See [IT[ §2, pp. 
57-58], 

We are interested in the prolongation of the solution germ to the right 
as long as possible in the domain T>. Actually the prolongation is unique, 
since the right-hand sides of ([T]) are analytic functions of r, m, P in T>. See 
P Chap.l, Sec.5], 

Especially we want to have a sufficient condition for that the prolongation 
turns out to be ‘monotone-short’ in the following sense: 

Definition 1 A solution (m(r),P(r)), 0 < r < r + , of (QJ ) is said to be 
monotone-short if r + < oo, dP/dr < 0 for 0 < r < r + and P —» 0 as 
r —)■ r + and if k+ > 0 and Q + > 0, where 


k + := lim K(r,m(r)) = 1 



2 Gm + A 


r —0 


Q + := lim Q(r,m(r), P(r)) — Gm 


with 


m + := lim m(r). 


r^-r+—Q 


Since the solution germ behaves as (l5i) as r —> +0, we assume 


( 6 ) 


A < — 5 — (p c + 3 P c /c 2 ) 
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in order that dP/dr < 0 at least for 0 < r C 1. 

Remark. If the condition (JHD does not hold but the equality 

(7) A = + 3 PJc 2 ) 

cr 

holds exactly, then the solution of (jT]) , (J5]) turns out to be (m(r),P(r)) = 

(47rp c r 3 /3, P c ), 0 < r < W3/L, where L := —— p c + A. In this very special 

c 2 

case, the metric is reduced to 

ds 2 = c 2 dt 2 — ^1 — ^r 2 j dr 2 — r 2 (d9 2 + sin 2 9d(f) 2 ) 

after a suitable change of the scale of t. The space at t = Const, is isometric 
to the half of the compact 3-dimensional hypersphere with radius \f3/L em¬ 
bedded in the 4-dimensional Euclidean space M 4 = {(^ 1 , £ 2 , £ 3 , £ 4 )} through 


£ 4 = r sin 9 cos 0, £ 2 = r sin 9 sin 0, £ 3 = rcos 6 ) , £ 4 



— r 


2 


Thus the horizon r = y/3/L — 0 is merely apparent. Of course the condition 
© is highly unstable. This is nothing but the ‘Einstein’s steady state inverse 
(1917)’ proposed in [7]. On the other hand, suppose that the inequality 

( 8 ) A > ^f(p c + 3 PJc 2 ) 

holds. Then the solution germ satisfies dP/dr > 0 for 0 < r -C 1. However 
it is possible that dP/dr become negative when prolonged to the right. This 
fact will be shown later. □ 


We introduce the variable u by 



dP 

p + P/c 2 ' 


Then we see 

u = —p 7 _ 1 0 . u (Ap 7 _ 1 /c 2 ), 

7-1 

p = Aiu~Vt p {u/ c 2 ), 

P = AAlu~flp(u/c 2 ), 
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where Q u ,Q p ,Q P are analytic functions on a neighborhood of [ 0 ,+oo[ such 

1 „ 1 

'7 - 1 \ — 


that f 2 u ( 0 ) = flp(O) = Q P (0) = 1 and A [ : = 


7 A 


The functions 


£l u , fL, Qp are depending upon only 7 and the function fh In fact we take 


n «(0 = > 


* n(o + ^C'Dn(C) 

1 1 + W) C ’ 
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fip(h) = ^(C)^u(O 7 " 1 with C = 


-770^(77) 


7 

Let us £x a small positive number <5^ such that these functions are defined 
and analytic on a neighborhood of [— +00 [. We put 

f Pc dP 7A 


u r 


' 0 P + p /c 2 7 - 1 


pZ-^ApV/c 2 ). 


2 Main result 


We claim 

Theorem 1 Suppose that 6/5 < 7 < 2. Then there exists a positive number 
e 0 (< 1 ) depending upon only 7 and the function such that if 

(9) u c < c 2 e 0 , A < (u c ) — e 0 , 

then the prolongation of the solution germ with |2|) to the right turns out to 
be monotone-short. 

This can be considered as the de Sitter version of the result by A. D. 
Rendall and G. B. Schmidt, [L5j. 


Let us sketch the proof. Using the variable u, we can write the Tolman- 
Oppenheimer-Volkoff-de Sitter equation (DO) as 

(10) c yy = Anr 2 Aifuifj^VLpiu/ c 2 ), 

du G { m +- y^r 3 

dr 9 U 2 Gm A 2 \ 

r v 1 - IFr ~H r ) 
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Here (w#) stands for max{«, 0}. Since we are assuming that 1 < 7 < 2 , which 
1 7 

implies // := - > 1 ,- — n + 1 > 2 , we can see that the functions 

7 — 1 7 — 1 

u n- (n # )^ +1 are of class C Yl (M). Keeping in mind it, we consider that 

the domain of the equation (jT 0 |) is 

D u = {(r, m, w)|0 < r, \m\ < 00 , — <5q < u/c 2 < + 00 , k > 0}. 

Let us perform the homologous transformation of the variables 

r = aR, m = a 3 b~ ■ AirAiM, u = bU , 

where a, b are positive parameters. We take b = u c and a which satisfies 

4:7rGAia 2 b~ = 1 . 


Let us write 

A := 


-A, 


AnGAi ’ 

Then the system (ITOT) turns out to 

dM 


a \= b/c 2 = u c /c 2 , /3 := b ~<~ l A = 


AnGAi 


(■ u c ) ■r-'A. 


( 11 ) 


dR 

dU 

dR 


= R 2 (U^n p (aU), 


R 2 


M + 1 — -aR 3 (U i )i^Q P (aU) - \$R} 
7_ ’1 

1 ~ 2a ii ~ l alSR2 ) 


Here ([/#) stands for max{[7, 0}. The domain of the system (TTT1) should be 
T>u = {(R, M , f/)|0 < R, \M\ < 00 , — 5q < U < 2, k > 0}, 
where, of course, 

M 1 9 

k = 1 — 2 a— — -a/3R 2 . 
ii 3 

Let us concentrate ourselves to0<a<l,0</3<l. We are considering a 
solution germ (M(R), U(R)) at R = +0 which satisfies 

r3 

(12) M(R) = n„(a) T + 0(R 5 ), 

u(R ) = 1 - (n„u,) + a n P (a) - flj A + o{r‘). 

We claim 
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Proposition 1 There is a positive number R 0 which depends upon only 7 
and the function 0 such that ( M(R),U(R )) exists and satisfies dU/dR < 0 
on 0 < R < Rq and (. M(R 0 ), U(Ro)) depends continuously on a, j3 G [0,1]. 

Proof is standard, and done by converting the system of differential equations 
(HU) to a system of integral equations under the condition (fT2j) as 


(13) 

g(R) 

U(R) 


"R 


R 3 
1 - 


O p (ct) 

R ±fi p (a)g(iO + 

1 - 2a\n p (ot)q(R')R' 2 - ±a/3i ?' 2 




Then, taking S sufficiently small uniformly on a, (3, we see that the mapping 
(q,U) ( q,U ), which is the right-hand side of the (TT3]) . is a contraction 

from $ = {(q, U ) G C[0, 5]|0 < q < C q , | < U <2} into itself with respect to 
a suitable functional distance, where 


C q := max 




< U < 2, 0 < a < 1 


For the details see HD pp. 57-58]. 


Let us come back to the proof of Theorem 1. The right-hand side of the 
system CD depends continuously on a, f3, and tends to (i? 2 ([/tt) 7-1 , — M/R 2 ) T 
as a —» 0, fi —> 0. The limit system 


dM 

dR 


R 2 (u*Y, 


dU_ 

dR 


is nothing but the Lane-Emden equation 


M 


(14) 


—A. 

~R?~dR 



Since we are assuming that 6/5 < 7 < 2 , say, 1 < ji < 5, the solution 
U = U(R) with U(0) = 1 of the Lane-Emden equation (fT4l) is short, that is, 
0 < U(R), dtJ/dR < 0 for 0 < R < & = ^( 7 ) and P( 6 ) = 0. See [5], 0. Of 
course we consider 







harmonically on R > . 

Thanks to Proposition 1, if e 0 is sufficiently small and if 


(15) 


a < e o and /3 < eo, 


then U = U(R ) exists and remains near to the orbit of U = U(R) on Rq < 
R < £i + 5r, 5r being small so that 


-y <U(^ + 5 r )<0. 


This is nothing but a direct application of j 6 j Theorem 7.4], Note that 
max{Qp(r/) | — Sq < q < 2 } depends upon only 7 and the function and 
we have — Sq < rj = all < 2 provided that a < e 0 < 1 and — 5q < U < 2. 
Especially if U(£i + 5r) < 0, then the radius R + of U(R) should be found 
in the interval ]0, + <5p[. This completes the proof of Theorem 1, since the 

condition (fT51) is nothing but flUJ). 

Note that (El) follows from (El) if e 0 is sufficiently small, since 



where 



is a function depending upon only 7 and the function and f 2 p+ 3 p/ c 2 ( 0 ) = 1 
so that we can assume that min{fi p+ 3 p/ c 2 ( 77 )| ~ 6^ < r] < 2} > e 0 . 

Remark. For the existence of u c satisfying (J9]) it is necessary that A 
enjoys 



Thus one may ask whether the real value of the cosmological constant of 
our universe satisfies it or not. But this question is not theoretical but 
experimental-observational and numerical. To answer to it is not a business 
of such a poor mathematician as the author of this article. □ 

Even if 7 < 6/5, the solution with the central density p c of the Tolman- 
Oppenheimer-Volkoff equation, that is, ([I]) with A = 0, or fl2J) , can be short, 
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if p c is large and the function P(p) is very much different from the exact 
7 -law for large p. For a sufficient condition for solutions to be short, see [121 
Proposition 3], the proof of [TT.> Theorem 1 ], Therefore we can consider such a 
case, supposing that 1 < 7 < 2 and the solution (m, P) = (m°(r), P°(r)), 0 < 
r < r°, of d2D with the same central density p c satisfies P°(r) —» 0 as r —> 
r []_ — 0, with being finite. Then the associated (m, u) = (m°(r), u°(r)), 0 < 
r < r° + , satisfies (TTU1) with A = 0, that is, 


(16) 


^7 = Anr 2 Ai(uf) 7-1 VL p (u/ c 2 ), 

du G ( m + ^ 3 AA7(M S )^rn P (n/c 2 )) 
dr /_ 2Gm\ 


r 2 11- 


c 2 r 


■) 


In order to extend (m°(r),u°(r)) onto r > r)J_, we put 
, 0 /„\ 


rrv(r) = rrv + {\= m u (r^)), 

Q/ c 2 (, /, 2 Gm° 


u J (r) = — ( log (1 


c 2 r o 


log (1 


2 Gm ( l 




; 2 r / 


for r > r°. Then the extended (m°(r),u°(r)) satishes (fT 6 jl on 0 < r < +cxd. 
Since the right-hand side of (ITU]) tends to that of (ITU]) as A —» 0 , the solution 
of (ITU]) under consideration exists and remains in a neighborhood of the 
orbit (m°(r), u°(r)) on 0 < r < r° + + 5 r , 5 r being a sufficiently small positive 
number, provided that A is sufficiently small. Thus we have 


Theorem 2 Suppose that the solution of the Tolman-Oppenheimer-Volkaff 
equation ([ID with central density p c is short. Then there exists a small 
positive number e\ such that, if A < e±, the solution germ of the Tolman- 
Oppenheimer-Voikoff-de Sitter equation m with the central density p c has a 
monotone-short prolongation. 

We should note that C\ may depend not only upon 7 and the function 
but also upon A, c, G and p c . In contrast with Theorem 1, we have no hope 
to specify the manner of dependence. 


3 Monotonicity 

Here let us give a remark on the monotonicity of the solution of the Tolman- 
Oppenheimer-Volkoff-de Sitter equation (P). 
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When we studied the Tolman-Oppenhcimer-Volkoff equation, that is, (CQ) 
with A = 0, or (J2]), we see that if ]0,r + [, r + < +oo, is the right maxi¬ 
mal interval of existence of the solution in the domain T>, then dP/dr < 0 
for 0 < r < r + and P — > 0 as r - > r+ — 0. Proof is given in ra, irn- 
In other words, we can say on the Tolman-Oppenheimer-Volkoff equation 
that, if the prolongation of the solution germ under consideration is short, 
it is necessarily monotone-short. However it is not the case on the Tolrnan- 
Oppenheimer-Volkoff-de Sitter equation with A > 0. Even if dP/dr < 0 for 
0 < r -C 1 under the assumption (J6]) , dP/dr may turn out to be positive 
during the prolongation. Let us show it. 

In order to fix the idea, we suppose u c — 1, and put 

c 2 

r = aR, m = a 3 ■ AttAiM, AttGAio 2 = 1, u = U, A = — A. 

AttGtAi 


Then the system (PQ) is reduced to 


(17) 


dM 

dR 

dU 


= R 2 u^ p (y/c 2 ), 

i 


dR + 


- wi? 3 )x 


X 


2 M 




c 2 R 3c 2 


-l 


Here n 1/(7 — 1). 

The right-hand side of the system (1171) depends continuously upon the 
speed of light c and tends to ^ R 2 U——(^M — — R 3 ^ as c —> 00 . This 
non-relativistic limit equation can be written as 


(18) 


1 d f p 2 dU \ 
H?dR\ dRJ 


U^-X. 


In this situation we are assuming that A depends upon c and c 2 A/(47rGAi) 
tends to A. Since (1T81) is the Lane-Emden equation when A = 0, we shall call 
it ‘the Lane-Emden-de Sitter equation’ supposing that A > 0. 

Although we are supposing 1 < /i < +oo(<^> 1 < 7 < 2), we observe the 
limiting case /i = 1(<=> 7 = 2). Then the equation (fl8|i is linear and the 
solution U = U(R) with U(0) = 1 is given by 


U(R) = A + (1 - A) 


sin R 
R 
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explicitly. We have 


U(R ) = 1 - -— -R 2 + 0{R a ) 

6 

as R ^ +0, and the condition (J 6 ]) reads A < 1. Suppose that \ < A < 1. 
Then we find that 

dU 1 — A / sin f? 

_ = _p o S fl__ 

turns out to be positive for 37 t/2 < R < 2tt and so on, while U(R) > 0 exists 
and oscillates on 0 < R < +oo, and converges to A as R —* +oo. Therefore, 
this explicit example tells us that, if 7 is near to 2 , c is sufficiently large, 
and c 2 A/(4:7tGAi ) is near to a number A in the interval [|, 1[, the behavior 
of the solution under consideration may be similar, that is, the prolongation 
of the solution germ with u c = 1 is not monotone. On the other hand, the 
condition (JSJ) reads A > 1. Then dU/dR > 0 for 0 < R 1 but dJJ/dR 
become negative and U(R) oscillates and tends to the limit A as R —> + 00 . 
Therefore, this explicit example tells us that, if 7 is near to 2, c is sufficiently 
large, and c 2 R/ (AnGAi) is near to A > 1, then P(r ) of the prolongation of 
the germ, which satisfies dP/dr > 0 as 0 < r <C 1, become decreasing. 

In the definition of ‘monotone-short’ solutions of the Tolman-Oppenheimer- 
Volkoff-de Sitter equation we have required that k + > 0 and Q + > 0. Let us 
spend few words concerning these conditions. 

Consider a solution ( m,P ) = (m(r),P(r)), 0 < r < r + , in T> such that 
dP/dr < 0, which requires Q(r,m(r), P(r)) > 0 for 0 < r < r + , and suppose 
P(r) —» 0 as r —y r + — 0, with r + being finite. 

When we are concerned with the Tolman-Oppenhcimer-Volkoff equation 
(J2J) with A = 0 , the condition k + > 0 follows automatically. Proof can be 
found in [ 12 j. Of course, if A = 0, then Q + = Gm + > 0 a priori. 

However if A > 0 it seems that we cannot exclude the possibility that 
k_i_ = 0 a priori. Generally speaking, since k > 0 in T), we have k+ > 0. 
Suppose k + = 0. Then 


2 Gm. 


A 






and, since k > 0 for r < r + and k + = 0, we see 


dn 


k j . : = 


2 G a 2Gm 1 

= Inn-5-47 rp H---- 

r=r + -0 r-Hr+- 0 c A & T z 


dr 
2 Gm+ 1 


-A r 


c 2 r 


2 - l Ar + = ^(!- Ar +) < °- 


On the other hand, we have 


Q + = Gm+ - 


c 2 A 


r 3 = £!r± (1 _ Ar 2 ) 


> 0 , 


since Q > 0 for r < r + . Thus it should be the case that 1 — Ari = 0 and 
Q + = 0. In other words, n + = 0 requires Q + = 0 and A = 1. This is very 
non-generic situation probably hard to occur, but at the moment we have no 
reason to exclude this possibility. 


4 Metric on the vacuum region 


Suppose that we have fixed a solution (m(r),P(r)), 0 < r < r + , of the 
Tolman-Oppenheimer-Volkoff-de Sitter equation fll]) which is monotone-short. 
Then we have the metric 


ds 2 = K + e 2u / c2 c 2 dt 2 - —dr 2 - r 2 du 2 

K 

on 0 < r < r + , where 

duj 2 = d6 2 + sin 2 ddcj) 2 . 

We should continue this metric to the exterior vacuum domain r >r + . Natu¬ 
rally, keeping in mind the Birkhoff theorem, we should take the Schwarzschild- 
de Sitter metric 


ds 2 


2Gm + 

c 2 r 



2 Gm + 
c 2 r 


^r 2 ') dr 2 — r 2 du 2 

O ' 


on r > r + . As the whole we take 


ds 2 = g 00 c 2 dt 2 — gndr 2 — r 2 du 2 , 
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where 




n + e 2 “( r )/ c2 (0 < r < r + ) 



-9 11 


(0 < r < r E ), 


with 


m(r) (0 < r < r + ) 
m+ (r+ < r < r^). 


Here the constants r E ,ri, (0 < 77 < r# < +oo), are the values of r of the 
so called ‘cosmological horizon’, ‘black hole horizon’, that is, k(v, m + ) > 0 if 
and only if 77 < r < r E - In other words, we have 


n(r,m + ) = — (r - r T )(r E - r)(r + 77 + r £ ). 


See [3j. But this situation is possible only if 


If (ED does not hold, then /c(r, m+) < 0 for all r > 0. However, since we are 
supposing k + = /c(r + ,m + ) > 0, the condition (TTUh is supposed to hold and 
we have rj < r + < r E . Let us discuss the regularity of this patched metric. 

First we observe the regularity of u(r). We claim 

Proposition 2 The function u(r ) is of class C 2 in a neighborhood, ofr + and 


u(r ) = B(r + — r)(l + 0(r + — r)) 


as r —* r + — 0, with B := Q + /r+K + . Hence p(r) is of class C 1 and 



Proof. Since u(r) satisfies the equation (fTUj) . whose right-hand side is a 
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that k + > 0. Therefore the continuous solution u(r) turns out to be of class 
C 2 and 

du _ Q+ _ 

dr r=r + -0 r+K + 

This completes the proof. 


Now we are going to see the regularity of g 0 o,gn- Since 


d 


dr 


-m\r) = 



(r < r + ) 

(r + < r < te) 

is of class C 1 , m(r ) is of class C 2 . Therefore gu is twice continuously differ¬ 
entiable across r = r + . 

Since u vanishes at r = r + — 0, goo is continuous thanks to the definition 
of k + . We see 


d 

~r9 oo 
dr 


2ki du 


r=r_|_—0 


c 2 dr 


r=r .|_—0 


2 Q+ 

c 2 r 2 


and 


T9 oo 
dr 


r=r ++0 


^2 Gm 


c 2 r 2 


2A 

_ T r 


r—r .|_ 


2Q+ 

c 2 ri' 


Therefore goo is continuously differentiable. We have 


dr 2 


9 oo 


4 n + / du \ '■ 


2k±( 


d 2 u\ 


r=r + -0 C 4 \dr J r=r +-0 C 2 \dr 2 Jr=r+- 0 


c 2 A , 2Q+ , 2(Q + ) 2 

I Q I 


But by a tedious calculation we have 
d 2 u 

dr 2 r=r + -o 

This can be derived by differentiating the right-hand side of the equation for 
du/dr, that is, the second equation of (TTUT) . Therefore we see 


r\u + 


9 4 2 

crriKi 


dr 2 


g oo 


r—r .|_—0 


d 2 

— TUdoo 

dr z 


r=r++0 


4 Q+ 

c 2 r 3 


2A. 


Hence goo is twice continuously differentiable across r = r + . Summing up, 
we have 


Theorem 3 Given a monotone-short solution of the Tolman-Oppenheimer- 
Volkoff-de Sitter equation we can extend the interior metric to the exte¬ 
rior Schwarzschild-de Sitter metric on the vacuum region with twice contin¬ 
uous differentiability. 
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5 Analytical property of the vacuum bound¬ 
ary 

Let us observe the analytical property of a monotone-short solution (m(r), P(r)), 
0 < r < r + , of the Tolman-Oppenhcimer-Volkoff-de Sitter equation ([TJ) . 
Proposition 2 tells us that the associated u(r ) belongs to C 2 ([0,r + ]) and 

u(r) = B(r + — r)(l + 0(r + — r )) 

as r —>■ r + — 0, where B = Q + /r+K + . Moreover we claim 

Theorem 4 Any monotone-short solution u(r), 0 < r < r+, of (QJj enjoys 
the behavior at r = r + — 0 such that 

"7 

u{r ) = P(r + — r)(l + [r + — r, (r + — i r)^ zrT ] 1 ), 


therefore 

pir) = ^ B ) J l ( r +-r)^ I (l + [r + -r, (r+- r)^]i). 

Pere [Ad, X 2 ]i stands for a convergent double power series of the form 

a klk 2 X^X k 2 \ 

fel+fe2>l 

1 7 

Proof. Let us denote p :=-so that - — n + 1. 

7 — 1 7 — 1 

First suppose that /i is an integer. Then proof is easy. In fact (m(r), u{r)) 
satisfies at least on 0 < r < r + the system of equations 

( 20 ) -—= Anr 2 Aiu^Vlpiu/c 2 ), 

dr 

du G^in + ^-r z AAfu^ +l Vlp{u/ c 2 ) j - ^r 3 

dr 9 2 Gm A 2 

rP 1--- r 

V c 2 r 3 

and (m(r),u(r)) —> (m + , 0) as r —> r + — 0. But, since // is supposed to be an 
integer, the right-hand side of the system (12U]I is analytic function of (r, m, u ) 
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in a neighborhood of (r + ,m + , 0). This guarantees that m(r),u(r) admit 
analytic prolongations beyond r = r + to the right, and completes the proof. 
Of course this analytic prolongation is different from the (7 2 -prolongation as 
a solution of (|I0|1 . since ^ 0) for u < 0. 

Now suppose that p is not an integer. Since u(r), 0 < r < r+, is monotone 
decreasing, it has the inverse function r = r{u ) dehned on 0 < u < u c such 
that r(u ) — > r + as u —* +0. Then we have a solution (m,r) = ( m(u), r(u )) 
of the system of equations 


(21a) 

(21b) 


dm 

du 

dr 

du 


—47t r 



2 Gm 
c 2 r 


1 - 


2 Gm 
c 2 r 



Q 1 ■ Aiu^Q^u/c 2 ), 



where 

Q = G [rn + ^r 3 AA7M M+1 0p(w/c 2 ) j ■ 

Since ( m(u),r(u )) —>• (m + ,r + ) and Q —> Q + > 0 as u —> +0 and the right- 
hand sides of fl21a|Hl21bj) are analytic functions of r on a neighbor¬ 

hood of (0,0, m + ,r + ), we can apply the following Lemma in order to get 


m(u) = m + + u[u, rt M ]o, 
r(u) = r + + u[u, u^]q. 

Here [-, -]o denotes a convergent double power series. 


Lemma 1 Let p > 1 and f a (x, x^, yi, y 2 ), a = 1,2, be analytic functions of 
x,x ll ,yi,y 2 on a neighborhood of (0, 0, 0, 0). Let (yi(x), y 2 (x)), 0 < x < 5, be 
the solution of the problem 

(22) = f a {x,x^ ) y 1 ,y 2 ), y a \x= 0 = 0, a = 1,2. 

Then there are analytic functions ip a of x, x M on a neighborhood of (0, 0) such 
that y a (x) = xip a (x, xJ‘) for 0 < x -C 1. 


A proof of this Lemma will be sketched in the Appendix. 
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Since dm/du ~ — Cu with C = A'Kr'/ r {n + /Q + )Ai, and du/dr —> —B, we 
have 


m = m A 


— Cu^ +l + 1Tii n U 
n> 2 


/m+1 


+ minW 

n>0,Z>2 




r = r+ ——u 
B 


n>l 


£ Q„y m+( . 

n>0,l>2 


If mi n 7^ 0 for 3n > 2, then dm/du would contain the term u^ n with n > 2. 
However it is impossible, since the right-hand side of (12lap cannot contain 
such a term. Therefore m\ n = 0 for Vn > 2. If Ci n 7^ 0 for 3n > 1, then dr/du 
would contain the term u^ n . However it is impossible, since the right-hand 
side of the equation (12 lbl) cannot contain such a term. Therefore C\ n = 0 for 
Vn > 1. Thus we have 


£ 


m = TH 4 - — Cu ^ +1 + > mi n u 

n>0,l>2 

1 


fin+l 


r = r+ ——u 
B 


n>0,l>2 


lin+l 


Moreover we can show that m/ n — cj n — 0 for 2 < / < n by induction on l. 
In fact, fix n > 2. Then m 2ri = C 2 n = 0, since, otherwise, dm/du, therefore 
the right-hand side of (121 all , or dr/du, therefore the right-hand side of (121b|) . 
would contain the term u^ n+1 with n > 2, which is impossible. Therefore 
m 2 n = C 2 n = 0. Consider 3 < l < n. Assume m n ' : i> = c n >j> = 0 for 2 < V < n' 
with n' < n,l' < l — 1. If m n i 7^ 0, or, c n i 7^ 0, then dm/du, or dr/du would 
contain the term u which is impossible by the induction assumption. 
Therefore m n i = c n i = 0 for 2 < l < n. This implies that 


r = r 1 


1 

B U 
1 


n>0,Z>n+l,Z>2 




(^+l)n+Z —n 


= r+ - —u( 1 + [u, u M+1 ]i). 

The inverse function u = u{r) then clearly enjoys an expansion of the form 
u = B(r + — r)(l + [r + — r, (r + — r) M+1 ]i). 


This completes the proof of Theorem 4. 
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Appendix 

Le us sketch a proof of Lemma 1. 

First we assume that p is not an integer but a rational number, say, 
/i = q/p,p,q G N, p > 2 and p,q are relatively prime. Note that a function 
given by a convergent power series 

= ^Cijx'ix 1 *) 3 , \cij\ < (6 < 1 ), 

can be rewritten as 

Q -1 /if 

<p(x) = X X c in x ' m+ \ \cin\ <-jf for 0 < n < p - 1, 

l n =0 

where 

cin = X(% | i + qJ = l,j = pJ + n, 3 J e N}, 

and M = M/(e5 p ),8 = <5/e. This rewriting is necessary, since not c^’s but 
C/jj's can be uniquely determined for the given function </?(x). 

In fact first we note that pn cannot be an integer for n = 1, • • • ,p — 1. 
(Proof: Let us deduce a contradiction supposing that nq/p is an integer. 
We can assume q < p, by, if necessary, replacing q by q' := q — [q/p\p. 
Since nq/p < n, we see that nq/p is either 1, • • • , or n — 1, therefore q/p is 
either 1/n, • • • , or (n — 1 )/n. Hence p is a divisor of n, a fortiori, p < n, a 
contradiction to n < p — 1, QED.) Hence pn + l = /in' + 1', n, n', 1,1' eN, 0 < 
n, n’ < p — 1, implies n = n',l = l'. Then we have a unique numbering 
{ n kJk)ke N °f ('«,/)’s such that /in^ + 4 < pn^+i + 4+i- By induction on k 
we can deduce = 0 for Vfc from ^ C[ n x im+l = 0 Vx. This means the 
uniqueness of the coefficients q n in the above expansion of </?(x). 

Anyway suppose 


p-i 


f a (x,x p ,y 1 ,y 2 ) = 

<3j lnk\k2 X -1 1 ^2 ) 


n =0 
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with 


® lnk\k 2 


< 


M 


^/+/Cl+/C2 


(0 < n < p — 1) 


and put 


p -1 


F(x, m , m) = E irarru E 

n=0 

M 1 - 1 


k 2 


Then the problem 

dY 

dx 

has a solution of the form 


1 — x/5 1 — 1 — yi/6 1 — ?/2/(^ 


= F(x,Y,Y), Y\ x=0 = 0 


Y = Mx( 1 + 
p-i 


M' 


EE 0 < Q n < — 


Z ra=0 




On the other hand (1221) has a formal power series solution 


p -1 


c « x tm+i 


y« = Y,Y, ( 

l n= 0 

where the coefficients c^’s are determined by a recursive formula 
— n r a = - h a 

'On 'L W+l,n / i i i ^Ini 

l + 1 + pn 


h a — ^ „a 
°LR — 


77“ r 1 . . . r 1 r 2 ■ ■ ■ r 2 

/ , a Znfcifc 2 0'(l)n'(l) c Z'(fci)ra'(fci) C Z"(l)n"(l) C l" (k 2 )n" (k 2 ) • 


Here the summation in the definition of 6“^ is taken over 


L = qJ + Z + Z'(l) + • • • + l\h) + Z"( 1) + • • • + l"{k 2 ) 

with /'(1), • • • , Z"(l), ■ ■ ■ > 1 and J G N, and 

pj + R = n + n'(l) +-b n'(ki) + n"{ 1) +-f- n // (/c 2 ). 

Then it can be shown inductively that |c^J < C/ n , which implies the conver¬ 
gence of the formal power series solution. This completes the proof. 

A proof by a similar and easier majorant argument can be done when /i 
is an irrational number. Let us omit the repetition. 
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